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We study models for superconductivity with two interactions: due to antiferromagnetic(AF) 
fluctuations and due to phonons, in a weak coupling approach to the high temperature supercon- 
ductivity. The nature of the two interactions are considerably different ; is positive and sharply 
peaked at (±7r,±7r) while is negative and peaked at (0, 0) due to weak phonon screening. We 
numerically find (a) weak BCS attraction is enough to have high critical temperature if a van Hove 
anomaly is at work, (b) V^^ (AF) is important to give d-wave superconductivity, (c) the gap order pa- 
rameter A(k) is constant(s-wave) at extremely overdope region and it changes to anisotropic s-wave 
as doping is reduced, (d) there exists a first order phase transition between d-wave and anisotropic 
s-wave gaps. These results are qualitatively in agreement with preceding works; they should be 
modified in the strongly underdope region by the presence of antiferromagnetic fluctuations and 
ensuing AF pseudogap. 
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One of the most significant works in recent condensed 
matter physics is undoubtedly the discovery of high Tc 

|. copper oxides [Q. It seems the high value of Tc can not 
be understood by the 'classical' BCS theory a large 

' amount of experimental and theoretical works have been 
performed. In fact, quite a few number of high Tc super- 
conductors (HTSC) are found. 

A phase diagram with doping and temperature is 
very rich, with the antiferromagnetic ordered phase, the 

[ superconductive phases - underdope and overdope-, a 
pseudogap region as well as a 'normal' metallic range. 
Electron doped compounds show a similar succession of 
AF, superconductive and normal metallic phases, with 
lower critical temperatures. It is notable that, in common 
with the classical BCS case, the superconductive state 
occurs through a condensation of cooper singlet pairing 

: (kT,-ki). 

Currently there are many works that ignore phonon- 
I mediated interactions and, in many cases, van Hove 
anomaly altogether. In these studies one considers 
mostly interaction due to AF fluctuations. The main 
justification for these approaches is that the energy scale 
for magnetic interaction is ~ lOOOfsT which is somewhat 
] larger than the Debye temperature of ~ 300K. In fact, 
the solutions of the gap equation with only positive in- 
' teraction due to AF fluctuations gives Tc ~ lOOK 

which is consistent with the experimental results. How- 
ever it is not surprising since Tc varies a lot depending 
on the coupling constants. 

On the other hand there are a number of works 
which show enhancement of Tc due to a van Hove 
anomaly. In most of these investigations, the main pur- 
pose is to obtain high critical temperature which is con- 
sistent with the Tc observed at optimal doping. Therefore 
almost all the works except ref . [6] assume that a van Hove 
anomaly is on the Fermi surface at optimal doping. In 



this way the decrease of Tc from the optimal as doping 
is increased or decreased is natural. However, the sys- 
tematic computation of doping dependence has not been 
done. A potentially serious problem of these approaches 
is that antiferromagnetic (AF) fluctuations which should 
be important in the underdoped region are not taken into 
account. 

Recently, Friedel and Kohmoto |0 consider the case 
where the van Hove anomaly is at 'Fermi level' for 
undoped (half-filling). Besides phonon mediated interac- 
tion, one from AF fluctuations is considered. Indeed, it 
was shown the AF instabilities decrease faster with in- 
creasing doping than the effect of the van Hove anomaly 
on superconductivity They argue the foUowings: 

Phonon or electron mediated weak BCS attraction is 
enough to have high critical temperature if a van Hove 
anomaly is at work near enough to the Fermi level. This 
could apply to electron doped compounds and also to 
compounds with Cu02 planes overdope in holes, where 
Tc decreases with increasing doping. If phonons domi- 
nate, it should lead to an anisotropic but mainly s su- 
perconductive gap and probably also in electron doped 
compounds. A d gap should develop as observed in a 
number of cases if AF fluctuations are strong enough. In 
the underdope range, the observed decrease of Tc with 
hole doping can be related in all cases to the develop- 
ment of antiferromagnetic fluctuations which produces a 
magnetic pseudogap, thus lowering the density of states 
at the Fermi level. The observed mainly d supercon- 
ductive gap then can be attributed to antiferromagnetic 
fluctuations. 

The essential feature of the BCS theory is the pair- 
ing picture. This has been well established in the 'clas- 
sical' superconductors. Also there is almost no sign of 
breakdown of this picture in HTSC. However, in or- 
der to obtain the so called 'BCS formula', one needs 
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the approximations: (a) The pairing interaction is 
weak, (b) The density of states is not too fast vary- 
ing near the Fermi surface, (c) The pairing interaction 
is constant within the cutoffs iksTo near the Fermi 
surface and zero otherwise. These approximations lead 
to the BCS formula: =i l.lSkBTo exp{-l/N{0)V}, 
A{T = 0)/kBTc ^ 1.75, the universal specific heat jump 
at Tc, the isotope effect with exponent a = 1/2 etc. 

We assume (a) throughout this work, but (b) is cer- 
tainly not appropriate as shown below. There is a large 
enhancement of the density of states close to the van Hove 
anomalies. The approximation (c) is oversimplified and 
will be replaced by weakly screened phonon mediated in- 
teraction, by AF repulsive interaction, and in addition 
by both. Therefore we expect that the BCS formula can 
be violated in experiments, even if the BCS theory with 
a singlet pairing is correct. 

In order to examine the above ideas, we consider a 
model on the square lattice with specific interactions be- 
tween electrons and solve the gap equation for general 
values of hole doping. The effects of pseudogap are not 
taken into account nor the correction due to the small 
couplings between Cu02 planes and the ensuing fluctua- 
tions of superconductivity near T^. Then Tc keeps rising 
while doping decreases. This contrasts with the behavior 
of HTSC whose Tc decreases and eventually vanishes in 
underdope region. 

The high Tc oxides have a square lattice of Cu atoms 
and O atoms in between the Cu atoms (here we neglect 
the orthorombic distortion which exists for some HTSC). 
We take an effective tight-binding picture to consider the 
Cu sites only, then with effective transfer t, we have the 
following results. 

For undoped compounds, the Fermi level is then a 
square. Doping by electrons and by holes produce al- 
most square Fermi surfaces with nearly symmetric devi- 
ations from the undoped square surface. This fundamen- 
tal symmetry between electron and hole dopings fits well 
the general symmetry of the phase diagram observed for 
electron and for hole doped compounds. 

As a result of this geometry the Fermi level sits also 
near a strong peak in the density of states, which di- 
verges logarithmic at (±7r,0), (0, ±7r) for the undoped 
compounds. This van Hove anomaly is characteristic of 
the (quasi) two-dimensional compounds, where it is much 
stronger than in more isotropic 3d compounds. 

The kinetic part of the square lattice model with 
only nearest-neighbor hopping can be written e{k) ~ 
—2t{coskx + cos ky) where the transfer t is typically 
0.25eV and energy will be measured in units of t from 
now on. In this energy dispersion, there are two saddle 
points which lead to the van Hove singularties at (±7r, 0) 
and (0, ±7r). The gap equation is 

A(k) = -iX:y..A(k')tanh(^)^, (1) 
where A(k) is the gap order parameter, Vkk' 



is the interaction b etween electrons, and Ek = 
y/ {e{k) — -|- A(k)2 where /U is the chemical poten- 
tial. 

-BCS isotropic interaction Vkk' = const. (= —3t) < 
. In the ovcrdopo range one can reasonably neglect the 
effects of AF fiuctuations in the 'normal' metallic range. 
A high critical temperature Tc can then be obtained in 
the mean field BCS approximation. 

The critical temperature is plotted in Fig. 1 for weak 
coupling region where N\V\ < 0.3 where N is DOS near 
the bottom of the band. At the bottom of the band 
Tc c:^ lOK which is quite reasonable for weak coupling 
BCS theory. It has an exact particle-hole duality. In the 
extremely overdope region, the BCS prediction is well 
satisfied because DOS is not so fast varying within the 
cutoffs. This is no longer true for lower dopings where 
DOS is larger and changes rapidly. The enhancement 
of Tc is considerable there, about 150K at undoped. The 
gap is s-wave since the interaction is constant. This result 
shows that weak BCS interaction is enough to have high 
T 

We stress that this enhancement of Tc is only induced 
by the DOS effect. The coupling does not play any role. 
It seems that it is not totally unrealistic to expect ex- 
tremely high Tc like room temperature. However, there 
exist AF effects near undoped in the all known high Tc 
superconductors. It is likely to suppresses Tc. 
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FIG. 1. Critical temperature Tc vs chemical potential for BCS 
attractive interaction (Vkk' = < 0)- 

- Weakly screened phonon-mediated interaction V^. 
Let us consider the interaction 

where q = k — k', gg is the electron-phonon coupling 

and 27r/qo is the screening length. This type of phonon 
contribution is used in ref.[8,9]. We take a large screen- 
ing radius about 30 lattice spacing because of poor 2d 
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screening. The sign of V^\/ is always negative. 

Critical temperature is displayed in Fig. 2. The gap 
A(k) now has an extended s-wave symmetry near half- 
filling. The new feature is that the gap not only depends 
on the angle of k but also on the absolute value |k| consid- 
erably being maximum at Fermi level. See Fig. 2 (a), (b), 
and (c) . Possibility of this phenomena has not been con- 
sidered seriously in the past. In most of the previous 
works on superconductivity, A(k) is computed for k at 
the Fermi level. 

A notable feature is that docs not decrease by dop- 
ing very much. In this region it is very hard to have pure 
HTSC samples. Thus it is likely that our results do not 
apply. For extremely doped region we have also high 
superconductivity. This behavior will be discussed else- 
where [ pr] . 
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FIG. 2. Critical temperature Tc for weakly screened 
phonon-mediated interaction (Vj^j^, = — |5(j|'^/(q^ +qo^)), and po- 
lar plot of gaps at T = for (a) /j = —3.6 (extended isotropic 
s-gap), (b) fi = —2.0 (slightly anisotropic s-gap), and (c) fi = —0.2 
(anisotropic s-gap). The distance from origin is the amplitude of 
gap for a given angle 9 = tan~'^{ky/kx). (d) (1)A{0) max / ksTc 
and {2)A{0)ave/kBTc vs /x. 

—AF fluctuation mediated interaction V'' . There are a 
number of attempt to obtain high Tc within the weak cou- 
pling BCS theory only from repulsive interactions which 
are originated from AF fluctuations Here we take 

the form of interaction by Monthoux et al. |^ This inter- 
action has peaks at (±7r, ±7r). We, however, take much 
weaker coupling about 1/10 compared with one used by 
Monthoux et al. [Q (The peak value of their interaction is 
about 80eV^). We obtain relatively low Tc {Tc 30K at 
undoped) and it vanishes at doping c ~ 0.1. The result 
obtained by using the coupling of ref. [4] is also confirmed. 

This contrasts with attractive interaction which always 
give superconductivity at low temperature for all doping. 
The gap is always d-wave for this type of repulsive inter- 
action 0]. Also we do not obtain a gap which breaks 
the square symmetry like the one studied by ref. [12] even 



in low temperatures. 

— The two types of interaction + V'^ . Finally we 
take both interactions into account. The two interactions 
do not cancel directly since has peaks at (±7r, ±7r) 
and at (0, 0). In any event, the peak of V^^., is fixed 
at (±7r, ±7r) and the effect of the interaction near Fermi 
level strongly depends on doping. On the other hand 
V^-^, does not so much depend on doping. Thus the rela- 
tive strength of the two interactions changes with doping. 

The results are shown in Fig. 3. In case both s and d 
gap can comprise the (mixed) solution of gap equation 
for a given point (/i, T) in phase space, the gap function 
is chosen as the one which has the lower free energy. A 
notable feature is that, although is not more effec- 
tive for superconductivity, it is very effective for d wave 
gap. We find a first order transition line between d and 
anisotropic s-wave gap. The latent heat is also displayed 
in Fig.3(d). There are experimental results which show 
a nodeless HTSC |l|,||. 
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FIG. 3. Critical temperatures T^and the first order phase tran- 
sition line (dotted line) for the interaction + V'' , and polar plot 
of gaps at T = for (a) /i = —3.6 (extended isotropic s-gap), 
(b) fj, = —1.0 (anisotropic s-gap), and (c) fj, = —0.2 (d-gap). 
(d) is the latent heat across the first order transition line, (e) 
{l)\A{0)\max/kBTc and {2)\A{0)\ave/kBTc vs ti. 

The quantity \A{0)\max/kBTc takes the usual BCS 
value about 1.75 for isotropic region (See Fig.3(el)). In 
the anisotropic region it is larger and about 2.5 which 
is not surprising because the gap is anisotropic and even 
has nodes in d-wave cases. This ratio with |A(0)|ai,e av- 
eraged over k is also plotted in Fig.3(e2). 

Our model has the particle-hole symmetry, thus we 
expect that the electron-doped compounds have similar 
gap structure. Mostly s-gap reported may be anisotropic 
and there may be a region of d-wave pairing in underdope 
region [p^ . 

We also find a novel scaling relation; Namely, the abso- 
lute value of average gap function |A(T)| as a function of 
T for different chemical potentials can be collapsed into 
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a universal curve by 

|A(T)| 
|A(0)| 



1 - 



= f{t/tcha 



where t = 1 — T/Tc is the reduced temperature and tchar 
is the characteristic reduced temperature for the decay 
of 1 — |A(T)|/|A(0)| as t increases for different chemical 
potentials, and / is a universal scaling function. Figure 
4 shows a good scaling collapse of the rescaled gap func- 
tions in terms of t /tchar for different chemical potentials. 
We observe that |A(T)| decreases as T increases by one 
fashion for /i below —1.4 and by the other fashion for fj, 
above —1.4 (See fig.4(a)). In fact Fig.4(b) shows that 
there are two different scaling functions for /j, below and 
above —1.4 at which Tc becomes minimum. This uni- 
versal scaling does not apply for the region of first order 
phase transition. 
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FIG. 4. The scaling collapse ofS-cfcaled gap functions versus 
t/ tchar foi' different chemical potentials from = —4.0 to ^ = 0.0. 
The inset (a) shows |A(T)| vs T for different chemical potentials, 
which show different decaying behavior as T below and above 
fi = —1.4. The curves (1),(2) for fi = -3.8,-3.6 decay in the 
same fashion, but differently from that of the curves (19), (20) for 
H = —0.2, 0.0. The inset (6) shows that there are two different 
scaling functions below and above /i = —1.4. 

To conclude, the presence of totally different types of 
interactions can be expected to lead to rich physics. 

We find that the phonon-mediated interaction is 
enough to have high temperature superconductivity. The 
role of AF fluctuation is to have the d-wave gap. More 
importantly it gives the AF pseudogap which we do not 
take into account by technical difficulties. Probably the 
most important prediction is the existence of the first 
order phase transition betwen d and anisotropic s gaps 
which may be tested in experiments with clean samples 
in the overdope region. 
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